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We analyzed time evolution of the localized charge in the system of two interacting single level 
quantum dots (QDs) coupled with the continuous spectrum states in the presence of electron-phonon 
interaction. 

We demonstrated that electron-phonon interaction leads to the increasing of localized charge 
relaxation rate. We also found that several time scales with different relaxation rates appear in 
the system in the case of non-resonant tunneling between the dots. We revealed the formation of 
oscillations in the filling numbers time evolution caused by the emission and adsorption processes 
of phonons. 
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I. INTRODUCTION 

Recent progress in the engineering and fabrication of 
well-defined artificial systems - quantum dots (QDs), 
leads to the possibility of ultra small electronic devices 
formation with a relatively high control of system pa- 
rameters (size, shape and energy spectrum) pj. In addi- 
tion to the QDs potential industrial applications, these 
nanoscale objects provide an ideal test ground for the 
study of basic physics including, many-body interaction 
effects, electron transport and time-dependent effects. 
Moreover, QDs integration in a little quantum circuits 
deals with careful analysis of non-equilibrium charge dis- 
tribution, relaxation processes and non-stationary ef- 
fects. These processes influence strongly on the electron 
transport through the system of QDs [1,11,13, H> II '0 • 
Electron transport in such systems is governed not only 
by the Coulomb interaction between localized electrons 
but also by the electron-phonon interaction 




Most of the theoretical works devoted to the problem of 
electron transport through the coupled QDs in the pres- 
ence of electron-phonon interaction deal with the tunnel- 
ing current and current-current correlations (shot noise) 
investigations [ll,[ll|]. Only a few attempts have been 
made to analyze phonon assisted localized charge relax- 
ation Qjl , El | , 15 1 . It was found theoretically that lateral 
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rates in parabolic QDs [12j |. Author considered the de- 
formation potential coupling between electrons and lon- 
gitudinal acoustic phonons, while neglecting piezoelectric 
coupling on the grounds of it's weaker contribution in two 
dimensional structures. On the other hand further theo- 
retical and experimental works suggested that electron- 
acoustic phonon scattering due to the piezoelectric field 
interaction is relevant for momentum and spin relaxation 
processes [l6l|.[l7| and may even provide charge decoher- 
ence in laterally coupled QDs [H,[l3]- In authors 
analyzed phonon induced single electron relaxation rates 
in the models of weakly confined single and vertically 
coupled QDs taking into account both mentioned above 
mechanisms. The regimes where each coupling mecha- 
nism prevails were found. 

In this paper we use the Keldysh diagram technique 
18] to analyze charge relaxation in a double QDs due 
to the coupling with the continuous spectrum states in 
the presence of electron-phonon interaction. Tunneling 
to the reservoir is possible only from one of the dots. 
We have found that electron-phonon interaction results 
in the increasing of localized charge relaxation rate and 
also leads to the formation of well resolved oscillations. 



II. THE SUGGESTED MODEL 

In the present paper we consider a system of coupled 
QDs with the single particle levels £\ and £2 connected 
with the continuous spectrum states. At the initial time 
two electrons with opposite spins are localized in the first 
QD on the energy level E\ (ni CT (0) = ran = 1). The second 
QD with the energy level £2 is connected with the reser- 
voir (e p ). Relaxation of the localized charge is governed 
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by the Hamiltonian: 

H = Hjj + H tun + H res . (1) 
The Hamiltonian Hp of interacting QDs 



In the case when initial charge is localized in the first 
QD and the second dot is empty, the third term in 
the eq.® can be neglected. Retarded Green's function 
Gf?(t',t) = [Gft T (*,*')]* yields density of states in the 
first QD and can be found exactly from the integral equa- 
tion: 



i=l,2a cr 

+ u; b + b + g(c+ a c 2a + c la 4 a )(b+ + b) (2) 

contains the spin-degenerate levels Si (indexes i = 1 
and i = 2 correspond to the first and to the sec- 
ond QD) and electron-phonon interaction. The cre- 
ation/annihilation of an electron with spin a — ±1 within 
the dot is denoted by c^/ciu and n a is the correspond- 
ing filling number operator. Operators b + jb describe the 
creation/annihilation of the phonons. ujq - is the optical 
phonon frequency and g - is the electron-phonon coupling 
constant. The interaction between the dots is described 
by the tunneling transfer amplitude T which is consid- 
ered to be independent of momentum and spin. 

The coupling between the second dot and the contin- 
uous spectrum states is described by the Hamiltonian: 



Hi, 



Cp a C 2o 



^pcr 



4J, 
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where t is the tunneling amplitude, which we assume to 
be independent on momentum and spin. By considering 
a constant density of states in the reservoir v , the tunnel 
rate 7 is defined as 7 = irvot 2 . 

The continuous spectrum states are modeled by the 
Hamiltonian: 



H r 



(4) 



pa 



where Cp a /cp a creates/ annihilates an electron with 
spin a and momentum p in the lead. We shall use 
Keldysh diagram technique to describe charge density 
relaxation processes in the considered system. Time evo- 
lution of the electron density in the QD is determined by 
the Keldysh Green function G^ which is connected with 
the localized state filling numbers in the following way: 



G^{t,t) = im{t) 



(5) 



Integro-differential equations for Green function 
Gy^(t,t ) without electron-phonon interaction has the 
form: 



G< T (M) 
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where Green's functions G\ R (t — t) and G 22 (t — t ) in 
the absence of coupling between the dots are determined 
by the expressions: 



G?f(*-*') = -*6(t - t ')e- fel <*-* ' } 

G22 (*-*') = - t ')e~ ie2(t " 4 ')-7(t-* ') ( 8 ) 

The eigenfrequencies Ei^ of equation ([7]) are deter- 
mined in the following way 

{E-e 1 ){E-e 2 + i 1 )-T 2 = Q 



£1,2 = + £ 2 - if) ± -£ 2 + *7) 2 + 4T 2 (9) 

Finally retarded Green's function can be written as: 



G R ?(t,i) = -iQ(t-t)( 
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and interaction with the continuous spectrum states is 
included in the Green's function G 22 (t — t ). Electron- 
phonon interaction results in the appearance of correc- 
tions to the Green's function G{\ T in the equations ([6]) 
and ([7]). Consequently the equation for Green function 
has the following form: 



r<R (+ +' \ — nOR . n^R^n^RnR _i_ nOR^R nR _i_ 
^llV 1 ! 1 ^ — Wi + Wi J ^22 Wi + ^nWi "+" 
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' 11 12 21 12 21 11 12 22 21 

(11) 

where self-energies Ef x (i,2; ), Ef 2 (i,t ), Y, 21 (t,t ) and 
E||(t, t ) can be written as: 



E^ 2 (i, t 
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In the cq. (fT2"|) the following ratio is considered 
(6) G$ 2 T {t ,t) = [G§ 2 T (t,t')}* and expression for Green's 
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function Gy£ '(t,t ) analogous to the equation (fT0|) has E® - are the eigenvalues of equation (fTS"]) with electron- 



the following form: 



phonon coupling constant 5 = 0. 



E 2~£l e -Ei{t-t) 

E\ — E2 
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) (13) 



The last three terms in eq. lfTTj) correspond to the next 



order perturbation theory in the parameter 



Con- 



sequently, localized charge relaxation in the presence of 
electron-phonon interaction is mostly governed by the 
term G?f Sf x Gf x . Acting with operator G°f _1 and con- 
sidering the terms with the accuracy ^2 one can re-write 
the equation in the following way: 



(G?f- 1 -T a G§| t -Ef 1 )Gf 1 (t J t') = 5{t-t) (14) 



The eigenvalues of equation ([14]) can be found from 
characteristic equation written in the following form: 



E 1.2 
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Equation for Keldysh Green function Gf 1 (i,t ), which 
determines localized charge time evolution ni(t) than has 
the form: 



G<(t,t) = G"<+G° 1 ?T 2 G^G< 1 + 

+Wl 2j ll < - r ll "+" Wl 2j ll < - r ll "+" Wl ^ll^ll 

where self-energy E^f^i ) can be written as: 



(20) 
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with phonon function D < (ti,t 2 ): 
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Consequently retarded Green's function Gf x can be 
written in the following form: 



Gf 1 {t,t) = ^-i@{t-t)A i 



-iEi(t-t ) 



(17) 
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acting with operator G^f 1 equation (l20l can be re- 
written in the following form: 



G° i n- 1 G? 1 (t,t) = (i? i -e 1 )Gf 1 (t,t) = 

POO 

= T 2 dt 1 -G° 22 H (t,t 1 )G< 1 (t 1 ,t') + 
Jo 

dt 1 [Z<(t,t 1 )G? 1 (h,t) + Sfi (i, ti)G& (ti , t')] (23) 



Green's function GJ^i, t ) = ini(t) is determined by 
the sum of homogeneous and inhomogeneous solutions: 



where - are eigenvalues of equation (fTS"]) . Coeffi- 
cients can be found from the system of linear equa- 
tions obtained in the first order perturbation theory in 
the parameter g 2 : 



5> 
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(24) 



Homogeneous solution of the equation can be written 
in the following way: 



(25) 



where coefficients A i correspond to the Green's func- 
tion Gyy, which is determined by the equation (|17p . 
Function G 22 \t\,t 2 ) can be written as: 



G< 2 T (h,t 2 )= E a ^ 

i'/=l,2 



e < • e 3 (26) 
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and coefficients •/ have the following form: 

iT 2 

an - a 22 = |£ o_ £ , f|2 

ai2 = = -an (27) 

Consequently one can find the inhomogeneous solution 
of the equation 

4 

A i- A r a i'j'- 

iji'j'=l 

-1 1 

, [e -i(Bj /+U0 )i _ e _ iBji][e i(B°;+"o)* _ eiB . t] (28) 
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Considering only the leading terms in parameters 

^j- in eq. (|28p , the following expression for the inhomoge- 
neous solution can be obtained: 

w o 7 ., ., 

j j =i 

.[e^Z+^'-e^*] (29) 

III. RESULTS AND DISCUSSION 

The behavior of filling numbers time evolution depends 
on the parameters of the system: energy levels detuning, 
the relation between tunneling rates and electron-phonon 
coupling constant, the value of optical phonon frequency. 
The general feature of all dependencies is the increasing 
of localized charge relaxation rate caused by the electron- 
phonon interaction. 

We start by discussing the filling numbers time evolu- 
tion in the case of the positive initial detuning between 
energy levels in the coupled QDs (Ae = E\ — e 2 > 0). Ob- 
tained calculation results are presented on the FigfT] It 
is clearly evident that electron-phonon interaction leads 
to the increasing of localized charge relaxation rate. The 
growth of the electron-phonon coupling constant g for a 
given set of system parameters results in the increasing 
of filling numbers relaxation rate. With the increasing 
of the initial detuning the influence of electron-phonon 
interaction on the charge time evolution is clearly pro- 
nounced (see FigJTkt). 

A critical value of the detuning exists in the system 
under investigation for a given set of parameters which 
corresponds to the relaxation regime changing. For the 
smaller values of the detuning charge relaxation takes 
place with the only one relaxation rate both in the pres- 
ence (see grey and black-dashed lines on the FigJTjb) and 
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FIG. 1: Fig.l Filling numbers time evolution in the presence 
of electron-phonon interaction in the case of positive initial 
detuning Ae. Black line corresponds to the case when g — 0, 
grey line describes the situation when g = 0.1 and black- 
dashed line - g = 0,2. a) Ae = 2,0, wo = 2,0; b)Ae = 1,0, 
ojo — 1,0. For all the figures values of parameters T = 0,6, 
7 = 1,0 are the same. 

in the absence of electron-phonon interaction (see black 
line on the FiglTJj) . The typical time scale which de- 
termines the localized charge relaxation is close to the 
value j res =2^—. For the larger values of detuning than 
the critical one, localized charge time evolution reveals 
two typical time intervals with different values of the re- 
laxation rates (see FigQJi). The first time interval re- 
laxation rate exceeds the relaxation rate of the second 
time interval both in the presence and in the absence 
of electron-phonon coupling. Without electron-phonon 
interaction the first time interval reveals charge relax- 
ation with the typical rate 7 res . The second time inter- 
val demonstrates charge time evolution with relaxation 

2 

rate close to 7„ orire s = Tres^jr- When electron phonon- 
coupling is involved, the filling numbers time evolution 
in the first time interval has the value 7 ~ 27 res and in 
the second time interval - 7 ~ 2 r y nonres . Consequently, 
electron-phonon interaction results in the two times in- 
creasing of localized charge relaxation rate. 

Let us now focus on the charge relaxation processes in 
the case of negative initial detuning between energy lev- 
els in the QDs (e\ < e 2 ) (see FigH]). One can clearly see 
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FIG. 2: Fig. 2 Filling numbers time evolution in the presence 
of electron-phonon interaction in the case of negative initial 
detuning Ae. Black-dashed line corresponds to the case when 
g — 0, grey line describes the situation when g = 0.1 and 
black line - g = 0, 2. a) Ae = -2,0, w = 2,0; b)Ae = 
— 1,0, ljo — 1,0. For all the figures values of parameters 
T = 0, 6, 7 = 1,0 are the same. The inset demonstrates 
localized charge relaxation in the case when g = 0,2. 

that in the case of negative detuning electron-phonon in- 
teraction also results in the increasing of localized charge 
relaxation rate, but this effect is less pronounced. Charge 
time evolution changes slightly in comparison with the 
situation when positive detuning occurs. In the case 
of negative detuning always exist several time intervals 
with different values of the relaxation rates (see Figf2]). 
For the small value of initial detuning (see FigJ^b) re- 
laxation rates on the both time intervals are very close 
to each other and to the value "/ res both in the pres- 
ence and in the absence of electron-phonon interaction. 



For the larger value of initial detuning in the case when 
electron-phonon coupling is absent the first time interval 
reveals charge relaxation with the typical rate j res and 
the second time interval corresponds to the relaxation 
rate jnonres- When electron-phonon coupling is consid- 
ered, relaxation rates increase slightly and they are con- 
tinue being very close o the values "f res and "fnonres for 
the first and second time intervals correspondingly. 

The most interesting effect in this energy levels con- 
figuration is the formation of oscillations in the filling 
numbers time evolution, caused by the emission and ad- 
sorption of phonons by the energy levels when detuning 
is close to the phonon frequency. Oscillations are well 
pronounced for the large value of initial detuning (see 
FigEK and the inset). For a given set of system parame- 
ters the oscillations amplitude is determined by the value 
of the electron-phonon coupling constant g (see Fig[2]). 
The presence of oscillations may even lead to the decreas- 
ing of filling numbers relaxation rate in comparison with 
the case when electron-phonon interaction is absent (see 
black-dashed and black lines on the FigfSk,) for a given 
set of system parameters in the particular time intervals. 



IV. CONCLUSION 

We investigated filling numbers time evolution in the 
system of two interacting QDs weakly coupled to the 
reservoir in the presence of electron-phonon interaction. 
It was shown that electron-phonon interaction results in 
the increasing of localized charge relaxation rate. The 
value of extantion is determined by the system parame- 
ters such as energy levels detuning, optical phonon fre- 
quency and the ratio between electron-phonon coupling 
constant and tunneling transfer amplitudes. 

We revealed that in the case of positive initial de- 
tuning between energy levels in the dots the influence 
of electron-phonon interaction is mostly pronounced and 
the relaxation rate increases with the growth of the initial 
detuning value. 

We found that when negative initial detuning is consid- 
ered, the influence of electron-phonon interaction leads to 
the formation of oscillations in the filling numbers time 
evolution. These oscillations are the result of phonons 
emission and adsorption between the energy levels. 

This work was supported by RFBR grants and by the 
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